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Fast NMR imaging was used to measure the density profile 
of a three-dimensional granular medium fluidized by vertical 
vibrations of the container. For container acceleration much 
larger than gravity, the rise in center of mass of the granular 
medium is found to scale as vq /N f3 with a = 1.0 ± 0.2 and 
13 — 0.5 ± 0.1, where vo is the vibration velocity and N is 
the number of grains in the container. This value for a is 
significantly less than found previously for experiments and 
simulations in one dimension (a = 2) and two dimensions 
(a = 1.3-1.5). 

PACS numbers: 45.70.-n, 81.05.Rm 



One technologically important way to fluidize a gran- 
ular medium like sand is to vibrate the walls containing 
the system. Basic issues include the scaling of the kinetic 
and potential energy of the grains with parameters such 
as the frequency and amplitude of the wall vibration, 
and the number of layers of grains. These scaling rela- 
tions are not yet fully understood for multidimensional 
systems, despite progress in simulations [jjj, experiments 
g, and theory |,|. 

We consider a granular system in gravity excited by pe- 
riodic vertical motion of the container z(t) = zocos(u>t). 
When the dimensionless acceleration T = zquj 2 / g is larger 
than one, the granular system enters a state of internal 
motion determined in part by T, the dimensionless fre- 
quency uj = uj{d/g) 1 / 2 , and the aspect ratio of the sys- 
tem. Here g is the acceleration of gravity and d is the 
grain diameter, so (d/g) 1 ^ 2 is the time for a grain to fall 
through its radius |l[ . For oJ <C 1 the container vibration 
induces large-scale structures in the granular medium 
such as density and surface waves. Surface structures 
are also observed for larger 57 « 1 in very low aspect-ratio 
systems that permit small- wavenumber symmetry break- 
ing m . Conversely, for ZU 3> 1 the container vibration is 
fast compared with the grain-scale motion. In this limit, 
the granular medium can assume a fluidized state analo- 
gous to a gas or liquid, and possibly amenable to theoret- 
ical methods developed for such systems The 
major distinguishing characteristic of a granular fluid is 
the continuous energy loss to inelastic collisions, which if 
sufficiently great can destroy the statistical uniformity of 
the system via inelastic collapse . 

In this paper we report an experimental study of a 
three-dimensional granular system in the fluidized regime 
ZJ w 4, r < 14. Past numerical and experimental 



work has concentrated on one- and two-dimensional sys- 
ljj ], due in part to limited computational 
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resources and the lack of experimental methods suit- 
able for examining the interiors of dense, flowing three- 
dimensional media. Recently, nuclear magnetic reso- 
nance (NMR) has emerged as a powerful noninvasive 
method for studying three-dimensional granular media 
|l5| [l7| . In the experiments reported here, NMR imag- 
ing was used to measure the density distribution during 
continuous vibration of the container. Incoherent grain 
motion limited image acquisition time to approximately 
1 ms, too short to permit ordinary two-dimensional imag- 
ing. Therefore one-dimensional images were acquired 
showing the density projected onto the vertical axis p(z). 
We find that p(z) deviates strongly from the exponential 
form e~ Cz expected for a dilute isothermal gas in gravity 
H . We focus here on the scaling of the gravitational po- 
tential energy, as measured by the center-of-mass height, 
with r and the number of grains N. In the present three- 
dimensional experiments the scaling is found to deviate 
even more strongly from ideal-gas predictions than was 
found in earlier two-dimensional studies |^|,^3) ■ 

The granular medium used for these experiments was 
composed of mustard seeds, of mean diameter 0.18 cm 
and mass 4.0 mg. Under magnification, the seeds ap- 
peared roughly spherical with typical eccentricity ±15%. 
Seeds are used because their oily centers give strong, 
long-lived NMR signals The seeds were held in 

a cylindrical container formed by a vertical glass tube 
of inside diameter 0.8 cm with a flat teflon bottom wall. 
The container was sufficiently tall to prevent collisions 
between the seeds and the top wall. It was determined 
visually that 18 ± 2 seeds were required to complete a 
monolayer in this container. The experiments were car- 
ried out under ambient atmospheric pressure. 

A vertical fiberglass tube was used to support the 
container at the center of the NMR probe, which was 
mounted into a vertical-bore superconducting solenoid. 
The static field was set to Bq — 1.00 T to reduce suscepti- 
bility contrast effects in these highly inhomogencous sam- 
ples. The lower end of the fiberglass tube was mounted to 
a loudspeaker, driven by a function generator and power 
amplifier, to provide vertical vibrations of the container. 
A flexible support near the top of the tube permitted 
the tube and container to move vertically while avoiding 
contact with the stationary NMR probe. 

The waveform and amplitude of vertical motion were 
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FIG. 1. Density profiles p{z) for N = 50 grains vibrated 
at 40 Hz. The acceleration of the vibration relative to grav- 
ity r was varied from to 14. The zero of the height axis 
z is arbitrary, as are the units of the density axis. With no 
vibration (r = 0), the grains are organized in distinct lay- 
ers which appear as peaks in p(z). At high values of V, the 
layer structure is eradicated and a smooth density profile is 
observed. 

continuously measured by a micromachined accelerome- 
ter (Analog Devices ADXL50) mounted to the fiberglass 
tube, close to the NMR probe. This accelerometer was 
initially calibrated against a more accurate unit, and its 
output was digitized and fitted. Data were taken at two 
values of the vibration frequency oj/2it — 50 Hz and 40 
Hz. Acceleration values were set to within ±0.05g. 

One-dimensional NMR images were obtained using a 
spin-echo sequence. Gradient pulses of fixed amplitude 
61.5 mT/m were used to encode vertical position infor- 
mation and an echo time TE = 1.0 ms was selected. 
Thus, the entire imaging sequence was complete within 
1.5 ms. These parameters were chosen to achieve ad- 
equate spatial resolution Sz = 800 /im within a time 
much smaller than the vibration period while holding sig- 
nal loss due to incoherent grain motion to an acceptably 
small level. This loss, measured by the echo size, was 
greater for smaller samples and larger vibration ampli- 
tudes. For a mid-size sample with N — 50 the loss was 
approximately 10% at the highest amplitude used. 

The NMR acquisition was triggered synchronously 
with the sample vibration. For each value of N and T, 
data for 100 different trigger delays spread uniformly over 
the vibration period were averaged together. Thus, the 
data presented here represent the sample density aver- 
aged over the vibration cycle. To remove effects of rf-ficld 



inhomogeneity the data were normalized to data obtained 
for a stationary water-filled tube. 

Figure |l| shows the vertical density profile as a func- 
tion of r for one of the (N, u>) points measured. With 
no vibration (r = 0) the profile is deeply corrugated, re- 
flecting the layering of grains in the stationary chamber. 
As the vibration amplitude is increased, the layer struc- 
ture is gradually smoothed out. This smearing is due to 
relative motion of the grains and not simply the vertical 
motion of the chamber, as the latter has an amplitude 
zq = (155 /im)r that is typically less than the grain di- 
ameter. 

Even at the highest vibration amplitude, it is clear 
from Fig. |l| that the density profile deviates strongly 
from the exponential form predicted for an isothermal 
ideal gas. Although the density tail at the upper (free) 
surface might reasonably be fit to an exponential, the 
density levels off and then decreases as the lower sur- 
face is approached. This agrees qualitatively with fea- 
tures observed in experiments and simulations for two- 
dimensional systems [HQ. 

To quantify the gravitational potential energy of the 
system, we have computed the center-of-mass height 
h cm = {f z P( z )dz)/(f p(z)dz). For each N, the limits 
of integration were kept fixed as T was varied to avoid 
introducing a systematic bias. Then the rise in center- 
of-mass height Ahem was computed by subtracting the 
r = value of h cm for the same N. Figure |2| shows 
Ah cm (T,N) for 50 Hz vibration frequency. The time- 
average height of the container is independent of vibra- 
tion amplitude, so positive Ah cm indicates an upwards 
displacement of the grains relative to the average con- 
tainer position. 

Small systematic dips can be noted in Fig. for specific 
values of N, for example N = 32 and N — 100. These do 
not reflect special states of the fluidized (r>l) state for 
these N values, but rather variations in the static (r = 0) 
packing of the particles which affect the subtraction used 
to compute Ah cm . These dips could be removed from 
Figs. |^ and ^ by assuming a power-law form for Ah cm 
and fitting to determine the T = limit of the fluidized 
state. We prefer to display the data after subtracting the 
measured T — h cm to avoid prejudging the functional 
form of Ah cm (T,N). 

To ascertain that functional form, we have plotted the 
Ah cm data scaled by various powers of the drive am- 
plitude zq and frequency ui. Figure |^ shows the best 
data collapse obtained in this way, which is found to be 
Ah C m/vo- Here vq — zquj is the velocity amplitude of 
the vibration. In this plot, the data for low acceleration 
r < 6 do not collapse while the data for high accelera- 
tion r > 6 all collapse onto a single function of N. The 
collapse is significantly better for Ah cm /va than it is for 
Ahcm/vy' 5 or Ah cm /vQ 5 , from which we conclude that 
Ah C m scales as Vq with a — 1.0 ± 0.2. As shown by the 
dashed line on Fig. 0, the variation of Ah cm with N at 
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FIG. 2. Rise in the center of mass of the grains relative 
to the unvibrated system, as a function of number of grains 
N and acceleration T. The legend gives the F value for each 
set of symbols. This data was taken with 50 Hz vibration 
frequency. 
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FIG. 3. The data of Fig. 2, divided by the vibration 
velocity vo and plotted on log axes. The symbols are the 
same as in Fig. 2. The eight curves for V — 7-14 lie nearly 
on top of one another, implying that Ah crn scales as Vq with 
awl. The dashed line is the function (0.055 s)/iV 1/2 . 



high r is approximately N 1 / 2 . From plots like this for 
a range of exponents we ascertain that Ah cm /vo oc 
with /3 = 0.5 ±0.1. 

A data set at fixed vibration frequency as shown in 
Fig. H can be used to determine the scaling with zq but 
not with to. When the data at vibration frequency 40 
Hz are added to the scaling plot it appears that velocity 
z uj is a better scaling variable than is the amplitude zq 
or the acceleration zquj 2 . However, due to the limited 
variation of u> this conclusion is not firmly established by 
our experiments. 

The simplest model for a vibro-fluidized granular 
medium is the elastic hard-sphere gas. In this model the 
mean grain velocity is proportional to the container vi- 
bration velocity, and the potential energy is proportional 
to the kinetic energy hence Ah C7n oc v^. This velocity 
scaling has been observed in experiments and simula- 
tions of one- dimensional granular systems (columns of 
nearly elastic beads) Two-dimensional experiments 

and simulations ]13| found Ah cm oc with signifi- 
cantly smaller a = 1.3-1.5, a result that is only partially 
understood jl|,0|. From a phenomenological viewpoint, 
our result a ~ 1 in three dimensions appears to continue 
the trend with increasing dimensionality away from elas- 
tic kinetic theory predictions. 

For a less superficial understanding we must consider 
the reasons for deviations from ideal-gas theory that have 



been proposed. The steady state of the granular system is 
determined by a balance between work done on the grains 
by collisions with the vibrating container bottom, and 
energy lost due to inelastic collisions between the grains 
0. (Other loss mechanisms discussed in the literature 
but not considered here include collisions with the side 
walls and viscous damping due to interstitial gas Q.) 

We first discuss the power Pi that is fed into the trans- 
lational degrees of freedom by collisions with the vibrat- 
ing bottom wall. An important dimcnsionless parameter 
is the ratio of the mean random grain velocity v g to the 
characteristic wall vibration velocity vq . For nearly elas- 
tic systems (restitution coefficient close to one) with a 
small number of particles per unit wall area, it is pos- 
sible to have v g /vo > 1. Conversely, less elastic sys- 
tems (including our experimental system) generally have 
v g/ v o *C 1 in the uniformly fluidized state at large U. 

Collisions between the grains and the vibrating bottom 
wall always transfer upward momentum to the grains, 
but only do positive work on the grains when the wall is 
moving upward. When v g /vo <C 1 the grains are unable 
to catch up with the wall on its downward stroke. In this 
case, equating the momentum transfer per unit time with 
the sample weight oc N gives the rigorous result Pj oc 
vqN Jl8| . The same result obtains for an asymmetric 
sawtooth vibration waveform (infinitely fast downward 
stroke) for any value of v 9 /vq. Thus the symmetry of 
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the drive waveform is only significant when v g /vo > 1 
For a symmetric waveform, it has been suggested 
and verified for two-dimensional simulations that Pi oc 
VqN f(v g /vo), where /(0) = 1 and f(x) is a decreasing 
function of x [Q. McNamara and Luding JlJ fit their 
simulation data for symmetric drive waveforms to f(x) = 
e~ Cx . We suggest their data may also be approximately 
represented by f{x) = 1/(1 + x). For large v g /v n , this 
gives Pi oc v 2 N/v g , a result that has been derived directly 
from kinetic theory for symmetric drive waveforms in the 
nearly elastic limit ||] . 

For a dilute, nearly elastic, isothermal system kinetic 
theory predicts that the density decays exponentially 
with height and the power lost to collisions P c oc N 2 v g 
Here isothermal denotes uniform granular tempera- 
ture T oc v 2 . In this limit potential energy is strictly pro- 
portional to kinetic energy, Ah cm oc T. Setting P c = Pi, 
for Vg/vo <C 1 or a sawtooth waveform the granular tem- 
perature scales as T oc v 2 , /N 2 . For symmetric waveforms 
with v g /vQ > 1 the scaling is T oc v 2 /N f|. The latter 
corresponds well to experiments and simulations for one- 
dimensional systems ]l2| , but not for two-dimensional 
systems p|,|l3| or the three-dimensional experiment re- 
ported here. 

Huntley Q] has suggested that correlations between 
grain motions at high density reduce the collision fre- 
quency below the kinetic-theory result for an ideal gas. 
The power lost to collisions is found to scale as P c oc 
N 3 / 2 v 2 . Combining this with the v g /vo ^> 1 result for Pi, 

Huntley found T oc v^ 3 /TV 1 / 3 . This is much closer to the 
experimental and simulation results for two-dimensional 
systems than is the unmodified kinetic theory prediction. 
If Huntley's expression for P c is set equal to Pi for the 
v g/ v o 1 condition appropriate to the present experi- 
ments, we obtain T oc vo/N 1 / 2 . This agrees well with 
the scaling seen in our experiments, suggesting that the 
degree of inelasticity may be the controlling factor for 
scaling exponents rather than the system dimension. 

Recently, McNamara and Luding JjJ carried out a se- 
ries of two-dimensional simulations in which Pi, P c , T, 
and Ah cm were measured as the relevant dimension- 
less parameters were varied. As expected, they found 
Ah cm oc T oc Vq in the limit of very low density that oc- 
curs at high drive velocity vq. For lower vq, they found a 
crossover to T oc v 3 / 2 which they traced to a reduction of 
P c below its kinetic-theory value when Ah cm is small. As 
the detailed mechanism for this reduction is unknown, it 
is not clear how these results can be applied to the present 
three-dimensional experiments. McNamara and Luding 
also found that Ah cm oc T is not well obeyed at high 
densities. 

The trend found in Ref. 0] of reduced power loss to 
collisions at high density agrees with Huntley's sugges- 
tion, although the functional form may be different. As 
the power input Pi appears reasonbly well understood, 



we also interpret our experiments as indicating reduced 
power loss to collisions P c for dense granular systems in 
three dimensions, as compared with ideal-gas kinetic the- 
ory predictions. 
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